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Abstract 

The large N infinite range spin glass is considered, in particular the number 
of spin components k needed to form the ground state and the sample-to- 
sample fluctuations in the Lagrange multiplier field on each site. The physical 
significance of k for the correlation functions is discussed. The difference 
between the large N and spherical spin glass is emphasized; a slight difference 

'^ . between the average Lagrange multiplier of the large N and spherical spin 

Q , glasses is derived, leading to a slight increase in the energy of the ground 

O I state compared to the naive expectation. Further, there is a change in the 

low energy density of excitations in the large A^ system. A form of level 
repulsion, similar to that found in random matrix theory, is found to exist in 

\^ • this system, surviving interactions. Even though the system is an interacting 

OO , one, a supersymmetric formalism is developed to deal with the problem of 

vQ I averaging over disorder. 

o 

0\ . I. INTRODUCTION 

-)— > 
<^ 

Spin glasses have been extensively investigated for a long time [0,0. Spin glasses with 

finite A^ (the definition of A^ is discussed below) are dealt with using a replica technique 

^ I which works wonderfully in the infinite range model. Unfortunately, it is known from other 

O ■ problems in random matrix theory that ^ that replica techniques sometimes fail, especially 

for nonperturbative properties. In this paper we look at the large A^ spin glass, and show 

that there exists a supersymmetric technique for dealing with this system. This system is a 

rS I disordered, interacting (although interacting only in the sense of the large A^ limit) system, 

c^ I that can still be dealt with in a fairly rigorous fashion. 

The large A^ generalization of the spin glass that we use is not new, although we will show 

that the behavior is much more complicated than has been found using replica techniques 

on this problem [^. For example, we define a number k, the number of spin components 

used to form the ground state, and find a non-trivial scaling of k with the system size; in 

the replica technique k is not even defined. We will find significant differences between the 

large A^ spin glass and finite A^ spin glass, much greater than the differences between large 

A^ and finite A^ unfrustrated systems. Many of these differences stem from the fact that 
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while an unfrustrated system will have the same ground state for large A^ as for finite A^, 
the frustrated large A^ systems will use a large number, k, of the available spin components; 
the finite N systems have far fewer spin components available. 

Of course, all physical glassy systems will have A^ = 1, 2, or 3, and will not be large A^ 
systems. The different behavior of the large A^ system in the glassy case makes the large A^ 
approximation less useful for glassy systems than it is for ordered systems. However, there 
may exist some optimization problems, or problems from computer science or combinatorics, 
which will lie in the large A^ limit, at least for finite V. Also, it may be possible someday 
to extend the techniques in this paper to the finite dimensional case, and obtain results on 
this system which cannot be obtained in the finite A^, finite dimensional problem. 

The infinite range model mentioned above is one of the simplest models of a spin glass. 
It is a system of V spins, all interacting with each either via Gaussian distributed random 
interactions. In the general case, we may let each spin have up to A^ components, so that 
the Hamiltonian is given by 

h = y: rA'jH,, (1) 

where z, j index the various sites and range from 1 to V , while /x indexes various components 
of the spin and ranges from 1 to A^. We impose the constraint that 

Y.(€? = N (2) 

We will take Hij to be a real symmetric matrix for most of the paper. In section XV we will 
discuss situations in which H may have complex entries and may also be complex. 

We will consider the problem in the large A^ limit. Consider the statistical mechanics of 
the system, integrating over all spin configurations with a weight e~^^ . It is well-known, 
and will be reviewed below, that in the large A^ limit the solution of a self-consistency 
equation will provide the properties of the system at any given inverse temperature (3 and 
for any given matrix Hij. A similar self-consistency equation will provide the ground state 
properties of the system. 

The problem that will be considered in this paper is the ground state and thermody- 
namic properties of the Hamiltonian given by equations ([1|,0), in the large A^ limit, when Hij 
is drawn from an ensemble of real, symmetric, Gaussian distributed matrices. There exist 
several solutions of such problems in the literature [^, using replica techniques. However, 
there are several problems with such solutions. Later in the paper we will discuss technical 
problems with that work. However, here, let us simply mention that two important prop- 
erties of the ground state (and the spin glass phase) are not addressed in that work. One 
property is the number of spin components used to form the ground state. Another is the 
strength of sample-to-sample fluctuations in the Lagrange multiplier field used to define the 
large A^ limit (discussed in a later section). 

Consider the number of components used to form the ground state. It is clear that 
a system with large A^ has more freedom to choose a ground state than a system with a 
smaller A^. There are more possible directions to point the spin components. To give a 
simple example, consider a system with three sites, and an anti-ferromagnetic interaction 
between the sites, so that Hij = 1 for alH, j. If A^ = 1, the lowest energy state is to take 



one site with spin +1 and two sites with spin —1, or vice- versa. If A^ = 2, one can find 
a lower energy state. For example, one can set 0i = -\/2(l,0), 02 = "\/2(— 1/2, v^/2), and 
03 = \/2(— 1/2, — a/3/2). The factor of \/2 in front of all the spins is simply to conform with 
the normalization chosen in equation (0); the important point is that the spins now use two 
different components. For the same system, with N > 2, one will still find that the spins 
only need to use two components. In fact, for this system, for any A^, the different spins will 
always span a 2 dimensional vector space. 

Let us define the number of spins used to form a given state. We will call this number 
k. We may define this number for any state as the dimension of the vector space spanned 
by the V different vectors 0f , where here as i = 1...V the spins 0f define V different vectors 
lying in the A^ dimensional space of vectors 0^. The vector space spanned by these vectors 
is a fc-dimensional subspace of this A^-dimensional space. 

The ground state of a system will have Goldstone modes, due to the possibility of rotating 
all the spins together. Ignoring these modes, a system will generically have a non-degenerate 
ground state; that is, all ground states will be linked by symmetry. In that case, the number 
of spins needed to form the ground state will be the dimension k of the space spanned by the 
vectors 0f used to form the ground state. However, we can imagine systems with degenerate 
ground states. For example, a system of V spins, with no interaction between the spins, has 
a very degenerate ground state manifold. All states have the same energy, so all states are 
ground states. The value of k for these states can be anything between 1 and V. When we 
refer to the number of spins needed to form the ground state, we will in this case use the 
smallest k of all the states on the ground state manifold. 

An important question is how big k is for the ground state in the large A^ limit, where 
the system can take k as large as it wants. One property we will find, due to the additional 
freedom to place spins in more directions for large A^, is the uniqueness of the solution to 
the thermodynamic self-consistency equation, and the absence of metastable states. This 
may be interpreted as an absence of replica symmetry breaking. 

The central result of the paper is the relation between V and the number of spin com- 
ponents k used to form the ground state. We will show that k{k + l)/2 < V for any given 
Hij. Then we will show that the infinite range system system, for large V, has k of order 
V'^^^. So, the bound is not saturated. 

Throughout the paper, we will be taking the large A^ limit first, and then considering 
systems with large, but finite V. By addressing the number of spin components used to 
given V, we will then be able to comment on the nature of the large A^ and large V limits. 
In fact, these limits are different, and the interchange of limits is not justified. One can 
guess that if A^ is of order V"^'^, or bigger, then the system will be in a regime described 
by large A^, while if A^ is smaller than V'^^^ the behavior will be completely different. It is 
likely that replica symmetry breaking will occur only if A^ is smaller than V'^^^. 

From results on random matrices 0, the idea of level repulsion is common. We will show 
in section II that for a given k, there will be a matrix which has k eigenvalues all of order 
1//3, so that in the large /3 limit these eigenvalues collide. This might seem to contradict 
the naive expectations of level repulsion; we will show in section V that a modified form of 
level repulsion still exists in the system. 

Another important result obtained in this paper is the magnitude of the sample-to-sample 



and site-to-site fluctuations in the Lagrange multiplier fleld Aj used to deflne the large A^ 
limit. This fleld also has a physical signiflcance, connected with the total energy of the bonds 
attached to a given site. This will provide information on the sample-to-sample variations 
in total energy. Within the replica formulations used previously, there are problems with 
calculating these fluctuations; this and other problems with the replica formulation will be 
discussed later. 

From a standpoint of formal technique, it is nice that we can use supersymmetric tech- 
niques. This is a glassy, interacting system, and it is useful to have such a system that 
can be treated without using replicas. Although the replica technique is very powerful, it 
sometimes has certain mathematical difficulties that a supersymmetric technique does not. 
For calculation of the spectra and correlators of random matrices, the Efetov supersymme- 
try techniques are superior 0. The supersymmetric technique used here is more closely 
connected to the supersymmetry introduced by Parisi and Sourlas [Q. However, we will 
make some connections to the Efetov supersymmetry technique in section XIII. 

The paper is set up as follows: first, we introduce the problem. In section II we review 
the solution of the large A^ problem via self-consistency. For any given if, we can exactly 
solve the problem by solving a self-consistency equation. The problem is to obtain results 
after averaging over different H . In section III and section IV we review theorems [^ showing 
that the solution of the self-consistency equation is unique and bounding the value of k. 

After this introduction and review of old results, we discuss some of what should be 
expected from the solution of the large A^ spin glass. We discuss level repulsion, paramagnetic 
and spin glass phases, and correlation functions. Section V introduces the idea of level 
repulsion in the interacting system. Section VI discusses the spherical model and section 
VII discusses the large A^ model in the paramagnetic phase. In the paramagnetic phase, 
these two models are equivalent in the thermodynamic limit. In the spin glass phase they are 
very different, and this difference is emphasized. In section VII we also discuss an analogy 
between k and the existence of Griffiths effects in finite dimensional systems. Section VIII 
discusses the consequences of the quantity k for the correlation functions and the physical 
significance of k. 

Then we proceed to the actual calculation for the large A^ system. In section IX the 
supersymmetric formalism is derived. In section X a simple example is discussed to help 
introduce some of the manipulations used in section XI when the supersymmetric formalism 
is used to obtain results on the ground state. In section XII we discuss the physical meaning 
of these results from section XI, and argue that the results of that section are also valid for 
the spin glass phase. 

Finally, we make some connections with Efetov supersymmetry and replica techniques, 
and discuss extensions of the problem. Section XIII shows some interesting connections 
between the supersymmetric formalism of this paper and the Efetov supersymmetry. Section 
XIV reviews problems with the replica solutions of the problem. Section XV discusses 
extensions of the results in this paper to more general matrices if, and section XVI is a 
conclusion. At the very end there are two appendices containing some of the mathematical 
manipulations needed for the calculations of section XL 

Finally, we will include a note on notation. Throughout, we will refer to vectors 0f or 
uf. The index i will range from 1...V and will index different sites in the problem. The 
index /z will range either from l...k or from 1...V and can be used to index different spin 



components on each site, or to index one oi k ot V different eigenvectors (these vectors he 
in the l^-dimensional space of vectors (pi). Eigenvalues will be written as E^. Lagrange 
multiplier fields which are a function of the site will be written as A,, Cj, and so on. 

II. LARGE N AND SELF-CONSISTENCY 

In this section, we will review the solution of the large A^ problem via a self-consistency 
equation. First we will discuss the solution of the thermodynamics of the problem at finite 
f3. Then we will discuss the self-consistency equation for the ground state. 

For the thermodynamics of the problem, we wish to compute the following partition 
function 



Z= J rf^fe-^^ (3) 

where 

h=y: m',H,, (4) 

Rewriting the constraint on the spins using a Lagrange multiplier field Aj, we can write the 
partition function as 

0frfAie-'^»(^^'^+^^^"^)'^^e^^' (5) 

where the integral over Aj extends from —ioo to +ioo. It is convenient to rescale Aj by a 
factor 13 to obtain 

I d4>f ciA,e-'^'(^(^'^+^^'^^^))'^^e^'^^^ (6) 

The next step is to integrate out the field 0, to derive an action for the Lagrange multiplier 
field. Then, after finding a saddle point for this field, it may be verified in the large N limit 
that fiuctuations about the saddle point are small. 

The partition function after integrating out is 

d\.e-'^^''^''(f^(^^^+^^^^^^'^e^^^' (7) 

The saddle point equation for Aj is then 

Gu = (3 (8) 

where we define the Green's function matrix G by 

{H,, + \)-' = G (9) 

In the large A^ limit, after solving for Aj, then the field can be treated as a free Gaussian 
field, with Hamiltonian 



f3{H,, + A,) (10) 

When solving for Aj, we must find a Aj such that 

Hij + Ai > (11) 

That is, H + Xi must be a positive definite matrix. 

The average energy of the system is given by the derivative of the log of the partition 
function with respect to /?. The result is 

V ^ 

-,-gA, (12) 

We can take a zero temperature (/? -^ oo) limit of the self-consistency equation (^ to 
find an equation defining the ground state of the system. Since equation (H) implies that 
Gii must diverge in this limit, some number of eigenvalues of if + Aj must go to zero. 

In the zero temperature limit of the problem, we can also solve the problem by giv- 
ing a configuration of spins, 0f which minimizes the energy. It is clear that if the given 
configuration is a local minimum of the energy, then for some Aj 

{H,, + XAjH^ = (13) 

for all fi. That is, the spin configuration must be made up of zero eigenvectors of Hij + XiSij. 
If the spins span a k dimensional space, then there must be at least k such zero eigenvectors. 
In fact, we can write the zero temperature problem as a problem of finding (pf and Aj 
such that equation ( p]3[ ) is satisfied, such that all spins have length equal to \/N, and such 
that 

Hij + XAj > (14) 

That is, the matrix Hij + XiSij must be positive semi-definite. 

To give a simple example in which k ^ 1, consider the system of three spins interacting 
anti-ferromagnetically with each other, as discussed in the introduction. We have Hij = 1 
for all i 7^ j. If we take Aj = 1 for all i, we find that H^j + Aj has two zero eigenvectors. 

We can choose to write these as \/|(l, -1/2, -1/2) and ^§(0, ^3/2, -a/3/2). By taking 



01 = V2(l,0), 02 = V2(-l/2, V3/2), and 03 = V2(-l/2, -^3/2), we have expressed the 
spins 0f as linear combinations of these eigenvectors, and satisfied the constraint on the 
length of the spin on each site. Here we have a system with k = 2. 

For large, but finite (3, the matrix H + Xi will have k eigenvalues which scale to zero as 
1//5. These are the eigenvalues which will vanish in the zero temperature limit. The other 
V — k eigenvalues will tend to non-zero limits in the limit of large f3. If these eigenvalues 
are denoted £^^ for fi = l...k, with eigenvectors ff, then an acceptable ground state is 

(15) 




In the thermodynamic hmit, it is possible for there to be a phase transition as the 
temperature is lowered. In this case, H -\- Xi will have k eigenvalues of order k/V. As a 
simple example, consider a system of V spins, such that Hij = —1/V for all i,j. This 
is a system with a ferromagnetic interaction between all the spins. For any /?, the self- 
consistency equation will be solved by taking Aj to be the same for all i. Then, H + Xi has 
V — 1 eigenvalues equal to A and one eigenvalue equal to A — 1. The self-consistency equation 
becomes 

^x-'+^{x-ir=p (16) 

For large V, this is solved in the high temperature, small (3, region by taking A = (3'^. As (3 
increases, however, there occurs a phase transition at (3 = 1. Past the phase transition, we 
must always have A > 1, in order for H + X to he positive definite. In fact, A will be equal 
to 1 + 0(1/1^). The self-consistency equation (|T6|) can be approximated by 

l + ^{X-ir' = f3 (17) 

We see that V — 1 eigenvalues remain gapped in the large V limit, while one eigenvalue 
becomes of order 1/V. Physically speaking, this is a state which is macroscopically occupied. 
For a frustrated system, in which k eigenvalues become macroscopically occupied, there will 
be k eigenvalues of order k/V. 

We expect that the properties of the system on the ordered side of the phase transition 
will be similar to the zero temperature properties. At zero temperature, k eigenvalues 
becomes equal to zero, vanishing as 1/(3. For finite, but large, (3, we have k eigenvalues 
of order k/V. For large V, the eigenvalues are all close to zero. So, it is reasonable to 
expect that the properties of the spin glass phase should be similar to the zero temperature 
problem. 

In the simple example of a ferromagnet, we had a gap between the one eigenvalue which 
was macroscopically occupied and the V — 1 eigenvalues above it. This will not remain true 
in the spin glass systems considered in this paper. There will be a spectrum of excitations 
above the k eigenvalues which are macroscopically occupied. This spectrum will have a 
small gap which vanishes as a power law in the large V limit; we will derive this gap in 
section XI and appendix A. All the states above the gap can be treated as a continuum and 
integrated over when obtaining the Green's function; there will be k states below the gap 
which must be treated more carefully. 



III. UNIQUENESS OF SELF-CONSISTENT SOLUTION 

Given the self-consistency equation defined in the previous section, we will show that the 
solution of this equation is unique. That is, we have the equation for all i 

Gu = 13 (18) 

where we define the Green's function matrix G by 



{Hij + A,)-^ = G (19) 

We will show that for any H there is only one solution Aj of this equation which has the 
property that Hij + Aj is a positive definite operator. 

The proof of this has been given before |Q, but we will review it again here. The proof 
proceeds as follows: assume for some given H, we have a A, that solves the self-consistency 
equation, with H + Xi positive definite. We will show below that for any infinitesimal 
change dH in H, we can find a change dA, in Aj such that the self-consistency equation 
is still obeyed. Also, we will show that the change in A, is well-behaved; that is, that the 
differential equation we will find leads to j^ always being finite. These two results imply 
that starting from any given H, we can deform H from the initial H to the point H = 0, 
and correspondingly deform Aj along this path, with the deformation of Aj along the path 
being unique and well-behaved. However, for H = 0, the only solution is clearly Ai = 1 
for all i. Now, if the initial H were to have two different solutions to the self-consistency 
equation with different Aj, then there would be a contradiction, since then we could find two 
different solutions to the self-consistency equation for H = hj deforming the two different 
solutions from the initial H to H = 0. So, the theorem will be proved. 

We will now demonstrate that -^^ always is finite, as stated above. Afterwards, we will 
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show that the deformations of H and Aj always keep H + Xi positive definite. 
Consider the equation 

Gu = 13 (20) 

Consider changes dH and dAj. Then we will have 

Gu + E G;MHik + dX,6,k)Gu = (3 (21) 

Then we have 

/ , GijdXjGji = — 2_^ GijdHjkGki (22) 

The right-hand side of this equation is some function of i. The left-hand side is some 
symmetric linear operator, acting on Xj. To prove the desired result we need to show that 
this linear operator is always invertible. This will follow if it can be shown that the linear 
operator is positive definite. To show that, we simply need to show that 

J2 dXiGijdX.Gj, > (23) 

for all vectors dAj. However, equation (^) is an immediate consequence of G being a positive 
definite operator. 

Now we will show that under these deformations G remains positive definite, which is 
equivalent to if + A remaining positive definite. It is impossible for G to acquire a zero 
eigenvalue since this would require that H + Xi have an infinite eigenvalue which cannot 
happen while H and A, remain finite. Then, the only way for G to acquire a negative 
eigenvalue is for one of the eigenvalues oi H + Xi to pass through zero. However, this would 
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require that one of the eigenvalues of G diverge. Looking at the self-consistency equation 
(§), and taking a trace, we find that 

Tr{G) = 13V (24) 



where V is the number of sites. Since all eigenvalues of G are positive, equation ( P^ implies 
that any given eigenvalue is bounded by V , so it is not possible for an eigenvalue to diverge. 
This concludes the proof. We have shown that, starting from a given H, and a solution 
of the self-consistency equation, it is always possible to deform to if = 0, while staying in 
the sector in which iiT + Aj is positive definite, and with the change in Aj always a regular 
function of the change in H. 

IV. MAXIMAL NUMBER OF COMPONENTS TO FORM THE GROUND STATE 

We will derive a bound on the maximal number of components needed to form the ground 
state of an arbitrary Hamiltonian for a system of V spins. The system will be in the large 
A^ limit so that it has as many components available as it needs. This theorem has been 
given before [^, but it will be reviewed here. We will consider here only the case of real 
Hamiltonians; others will be discussed in the last section. The result for real Hamiltonians 
is that fc, the number of components needed, satisfies the inequality k{k + l)/2 < V . At 
the end of this section we will explicitly construct a sequence of systems which saturate this 
bound for any k. 

The proof of the inequality k{k + l)/2 < V is very simple. Suppose the system needs 
k spin components to form the ground state. Then, for small perturbations 5H of the 
Hamiltonian H, the system will continue to need k spin components to form the ground state. 
However, as shown above in the discussion of the large A^ self-consistency equation, this 
means that the operator H + Xi has at least k zero eigenvalues. For a real, symmetric matrix, 
such as H + Aj, to have k eigenvalues requires adjusting a total of k{k + l)/2 parameters. 
Since we have assumed that for any small perturbation of H the system continues to have 
k zero eigenvalues, the only parameters available to tune to maintain the k zero eigenvalues 
are the V different values of Aj. Therefore, we have that k{k + l)/2 < V. 

Note that this result is not a result for an ensemble oi H; it is a result for any H. To 
express the proof more formally, let if be a V-hj-V matrix. The dimension of the space of 
6H is V(y + l)/2. The space of matrices {H + 6H + Aj + 6\i) that satisfy equation (g) has 
dimension V(y + l)/2 — V since we have V parameters 6\i to adjust to satisfy the V self- 
consistency equations. However, the space of V-hj-V matrices with k zero eigenvalues has 
dimension V{V + l)/2 — k{k + l)/2. If we need to use k spin components to form the ground 
state, then for small SH every matrix {H + 6H + Xi + 6Xi) has k zero eigenvalues. So a space 
of dimension V{V + l)/2 — \^ is a subspace of a space of dimension V{V -|- l)/2 — k{k + l)/2; 
this is only possible if k{k + l)/2 < V. 

For non-generic Hamiltonians H, we can have k zero eigenvalues of if + Aj, with k{k + 
l)/2 > V. However, this requires also adjusting some parameters in the matrix H to produce 
these extra eigenvalues. Since we have assumed above that for arbitrary small perturbations 
of H there will still be k eigenvalues, we have excluded this non-generic case. The case in 
which k{k + l)/2 > V^ is a case in which the Hamiltonian has a degenerate ground state, 
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as discussed in the introduction. However, the number of spins needed to form the ground 
state will still satisfy k{k + l)/2 < V. 

Finally, we can check the total number of parameters available to satisfy the self- 
consistency equation. To satisfy the self-consistency equation we have to satisfy a total 
of V different equations, one for each site. We have V parameters Aj available to satisfy the 
equations, so the number of parameters and equations balance. We have used k{k + l)/2 of 
these parameters to produce the required k zero eigenvalues, so there are V — k{k + l)/2 of 
these parameters left to satisfy equation (§). Given these k zero eigenvalues, consider the 
space of all vectors 0f , such that {H + Xi)ij(pj = 0. Here i = 1...V and fi = l...k. This 
space is of dimension k"^. There is a k{k — l)/2 dimensional group of rotations of this space; 
this leaves a k{k + l)/2 dimensional space of distinct states. Any state in this space is a 
good ground state, so long as it satisfies the constraint on the length of the spins. So, we 
have a total oi V — k{k + l)/2 parameters from the different Aj and a total of k{k + l)/2 
parameters from the possible 0, so there are V different parameters available to satisfy the 
self-consistency equation. 

To give a simple example of this balancing of parameters, consider the case k = 1. Then, 
one parameter is required to produce a single zero eigenvalue. There are V — 1 parameters 
left over, which are sufficient to make the eigenvector 0j have the same squared amplitude 
on each site. Then, we have one parameter available, which is the magnitude of (pi, to make 
this amphtude on each site equal to unity. 

One immediate consequence of this that for sufficiently large A^, the system cannot have 
any metastable states. Metastable states are spin configurations which are local minima, 
but not global minima. We have shown above that for any V, we can construct a ground 
state using k spin components where k{k + l)/2 < V. More generally, we have shown that 
any local minimum can be constructed using k spin components with k{k + l)/2 < V. Let 
us suppose we have an A^ component system, where A^ is equal to 2kmax where kmax is the 
largest k which obeys k{k + l)/2 < V. Suppose we have a local minimum of the system, 0f , 
which is distinct from the ground state pf. Since we have 2kmax total spin components, we 
can arrange things so that (pf is nonvanishing only for fi = l...kmaxi while pf is nonvanishing 
only for /i = k^ax + ^■■■2kmax- Then, consider the state yjl — 5^0 + 5 p. As 5 is changed 
from to 1, this provides a path starting from such that, near (5 = 0, the change in energy 
vanishes to first order in 5 but is negative to second order in 5. This implies that in fact 0j 
could not have been a local minimum, and that for A^ > 2kmax there are no local minima 
other than the ground state. 

Given the absence of metastable states shown in this section, and the uniqueness of 
the self-consistency equation shown in the last section, it is reasonable to suppose that 
replica symmetry breaking is absent in the system. Of course, replica symmetry breaking is 
only defined as a concept in the replica formalism, but the various concepts of inequivalent 
thermodynamic states, and so on, should be absent in the large A^ problem. This may be 
the reason that replica treatments of the large A^ problem do not find any instability about 
the replica symmetric solution ||^. 

Finally, we would like to show that the bound derived above for the maximum fc as a 
function of V is sharp. That is, for any k, we will construct a system using V = k{k + l)/2 
spins that employs k spin components to form the ground state. The procedure is discussed 
below and illustrated in figure (0). 
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First let us look at systems with small k. For A; = 1, we have k{k + l)/2 = 1. Obviously, 
a system consisting of one spin always employs one spin component. For A; = 2, we have 
k{k + l)/2 = 3. A system of three spins, all interacting anti-ferromagnetically with each 
other, with Hij = 1 for all i ^ j, will have a ground state using two spin components. This 
example was discussed in the introduction. 

To find systems with higher k, we will proceed inductively. Assume that, for some k, there 
is an Hij acting on a system of V^ = k{k + l)/2 spins such that the ground state uses k spin 
components. We will construct a Hamiltonian acting on a system of V^+i = {k + l)(A; + 2)/2 
spins which uses k + 1 spin components. Note that Vk+i — Vk = k + 1. That is, we have an 
additional k + 1 spins to use. 

Let the ground state of the system of Vk spins be given by pf where i = l...\4. Assume, 
without loss of generality, that the vectors pi for i = l...k span a /c-dimensional vector space. 
Then, consider the following Hamiltonian 

Vk k k 

H=Y: <t^f<P^H,, + Y: ct>tK,+^ + e E('/'f + 0y.+.)0v.+i (25) 

j,jr' = l i=l j = l 

where e is taken to be sufficiently small. 

The above Hamiltonian can be thought of as follows: there are still Vk spins with the 
same interaction as before. These spins are those with i = l...\4. An additional A; + 1 spins 
have been added. First we will consider the system defined by adding only the first k of 
these spins, then we will consider the system defined by adding all A; + 1 of these spins. 

We make each of the first k of these spins (these are the spins with i = Vk + l...Vk + k) 
have an anti-ferromagnetic interaction with one of the spins in the first Vk spins. Let us 
consider this system consisting only of the original Vk spins and these first k additional 
spins. This system has Vk + k = Vk+i — 1 spins. The ground state of this system is clear. 
The first Vk spins are in the ground state of the original system (0f = pf ), and the next k 
spins each point opposite to one of the spins in the first Vk spins. This produces k pairs of 
spins, such that the spins in a pair point in opposite directions. For example, 0i = —(py _,_i, 
02 = -0\4+2> and so on. 

Finally, we add the last spin, to produce a total of V^+i spins. This one additional spin 
has a weak anti-ferromagnetic interaction with the 2A; spins in these k pairs. 

Starting with the first Vk+k spins fixed in their ground state, adding the V^+i-th spin does 
not change the energy at all because the anti-ferromagnetic interactions cancel. However, if 
in a given pair of spins, one can slightly bend the pair, so that the two spins form an angle 
of less than vr, then one can gain energy from adding the V^+i-th spin. Since e is small, the 
energy gain is small. However, if the angle between two spins in a pair is n — 6, than the 
energy gain is first order in 6, while the energy cost is second order in 6. So, even for small 
e, it is advantageous to distort the pairs of spins. As chosen above, the first k spins span a 
A;-dimensional space. 

To distort a pair and gain energy, each spin in that pair has to have some component of 
its distortion in a direction opposite to the direction in which the V^+i-th spin lies. Also, for 
small distortions the direction in which a spin distorts must be orthogonal to the direction in 
which it points, since all spins must have fixed length. So, since each pair must distort in a 
direction orthogonal to that in which it lay orignally, and (pv^+i must have some component 
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in the direction in which that pair distorts, we find that the V^+i-th spin must have some 
component which is orthogonal to 0j for i = l...k. Then, the system of Vk+i spins must span 
a k + 1-dimensional space, as needed. 

In figure (0), we illustrate this procedure, in the simple case of the transition from k = 1 
to k = 2, going from l^ = 1 to V = 3. A solid arrow is used to illustrate the spin in the 
original system (the k = 1,V = 1 system), a dashed arrow is used to illustrate the spin 
which aligns opposite to that spin, a dotted arrow is used to illustrate the V^+i-th spin (the 
third spin), and small arrows are used to illustrate the direction of distortion of the first two 
spins after the third spin is added. 

V. LEVEL REPULSION 

An interesting question is the existence of level repulsion for the large A^ system. We 
know that if Hij is drawn from an ensemble of matrices that its eigenvalues will repel each 
other [^ . This repulsion is for a noninteracting system. This can be seen as the result of a 
Jacobian by writing 

J dHe-''^^'''" = J dEdOlllE,- E.llle-''^^''^" (26) 

i<j i=l 

where we have written H = O^ EO with O an orthogonal matrix and E a diagonal matrix 
of eigenvalues Ei of H . The factor \Ei — Ej\ causes the eigenvalues to repel each other. This 
can also be stated by saying that requiring two eigenvalues of H to be equal to each other 
requires tuning two parameters of if, instead of one parameter as might naively be thought. 
However, for the interacting large A^ system we are considering, the eigenvalues of if + Aj 
will be different from those of H. In fact, we know that the solution of the self-consistency 
equation can force k of those eigenvalues to scale to zero as 1/13. This means that the 
eigenvalues become very close to each other, differing only by an amount of order 1//5. Let 
us denote the eigenvalues of if + Aj by E^, where /x = 1...V . Let us have the first k of these 
eigenvalues be the ones that scale to zero. Among those k eigenvalues, there still is a kind of 
level repulsion. Consider a problem in the limit (3 —>■ oo. Then, the solution to the problem 
is a configuration of spins 0f defining the ground state of the system. These spins span a 
A;- dimensional space and so we may assume that /i = l...k. By a rotation of the spins, we 
can further assume that 0f is orthogonal to 0f for n ^ u. That is 

T.^t<t^'t = ^ (27) 

Then, by permutation of the spin components /i we can arrange it that 

T.m' = ^ (28) 

As /? — i> cx), S^ — i> for /i = l...k but I3E^ stays non-zero. 

Let us show the existence of level repulsion among the E^ in a simple model. Consider a 
system of three spins, interacting with some matrix Hij. We know that if Hij = 1 for all i y^ j 

12 



then we have the antiferromagnetic system discussed in the introduction. It can be checked 
that the particular ground state of this system given in the introduction satisfies equation 
(0) and that in this case, where 0,^ = ^2(1, -1/2, -1/2) and 0^ = ^2(0, ^3/2, -v^/2) 

h<Plf = h<f>lf (29) 

So then 

/3Ei = PE2 (30) 

This system has k = 2 and both of the eigenvalues which go to zero are equal. So, we have 
two levels overlapping. We would like to know how many parameters are needed to guarantee 
that the two levels overlap. If only one parameter is needed, then there is no level repulsion. 
If two parameters are needed, then there is level repulsion. The parameters we can vary 
are the different elements of Hij] there are three of these, since H is symmetric and the 
diagonal elements of H are unimportant. However, in fact there are only two independent 
parameters for determining the ground state of the system, since multiplying all elements 
of i7 by a common, positive factor does not change the ground state. 

Even if we do not rotate the spin components to make equation ( pTl) true, a simple 
equation for Ei, E^ can still be derived. It can be shown that 



1 1 
El E2 




(31) 



Consider matrices H near the uniform antiferromagnetic H considered above. There are 
two independent parameters defining the matrix H, if all we are interested in is ground state 
properties. There are also two independent parameters defining any set of states equivalent 
under rotations, since there are three spins and one overall rotation mode of the three 
spins. Let us choose the overall rotation so that we continue to pick ground states with 
01 = -\/2(l,0). Then it is easy to show that small variations in parameters in H produces 
small changes in 02, 03, with no singularities in the Jacobian relating changes in the two 
independent parameters of H to the two parameters defining 02, 03. So, instead of asking 
how many parameters of H must be tuned to produce Ei = E2, we can ask, when considering 
level repulsion, how many of the two independent parameters defining 02, 03 must be tuned 
to produce Ei = E2. Looking at equation (0), we see that we must have 

T.i^]f - E(<^')' = (32) 

i i 

and also 

E20.V? = O (33) 

i 

This requires tuning both available parameters to satisfy the two requirements. 

Then, level repulsion assume an interesting form. Let us forget about the underlying 
Hamiltonian and simply consider arbitrary spin configurations 0f . If we had forgotten about 

13 



the orthogonality requirement of equation (pT]) and simply used equation (pH]), it would only 
require one parameter to make two levels coincide. However, either if we rotate the spins to 
satisfy equation ( pTj) and then use equation ([28|), or if we directly use equation (pT]), we find 
that it requires additional parameters to make levels coincide, and this implies that there 
still is a form of level repulsion in the interacting system. 

VI. SPHERICAL MODEL 

We will discuss a simplification of the large A^ spin glass. This is the spherical model. This 
will be a useful simple example. Here the system will always condense into one eigenvalue 
(fc = 1). The spherical model is defined by using the Hamiltonian 

H = J2(j),<PjH,^ (34) 



«j 



subject to the constraint that 



Y.i<t>^r = V (35) 



Here, there is only one component of (f) on each site, but the constraint is relaxed to the 
requirement that the sum of spins over all sites be equal to V . The problem is again solved 
by self-consistency, with the self-consistency equation (again, we rescale A by a factor of (3) 

Tr(G) = (5V (36) 

where 

G={H + \)-^ (37) 

where A is now independent of the site index i. 

Consider an ensemble of matrices Hij given by the measure dif e"^"^""^^ , where r is 
some number of order unity. Then it may be shown that, averaged over different Hij, we 
have IHl 



Gii = 2rA - 2Vr2A2 - r (38) 

Also, we have that, on average. 



Tr(G) = 2VtX - 2VVt^X^ - r (39) 

Using this in equation (p6D, we find that 

'-hi («) 

One way of obtaining equation (^) is by knowing the spectrum of H. It is known 
that H has a distribution of eigenvalues obeying the Wigner semicircle, shown in figure (^. 
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The eigenvalues, in the large V limit, range from —y^ to JK The density of eigenvalues is 
given by Vp{E), where p{E) = —J- — E'^. This implies that 



TrG= / dEVp{E)^-^ (41) 



Doing this integral yields equation (|39D . 

At /5 = 2y^, we have A = J^ from equation (^). At this point the spectrum of if + A 
becomes gapless, and a phase transition occurs. The sum over eigenvalues needed to define 
TrG can no longer be approximated by equation (PT|), and now must separately include the 



contribution from the one lowest eigenvalue of H . So, we have that 



TtG= \ + / dEVp{E)-^ (42) 



where Emin is the lowest eigenvalue of H. The self-consistency equation (|^) can be satisfied 
in the spin glass phase only if A + Emin is of order 1/V. In this phase, we will have A ~ -, 
so that H + \ will be approximately gapless. 

More precisely, it is known 0,^ that the lowest eigenvalue of the system, Emm will lie 
within a tail that extends to a distance of order V~^'^ below -. Then, A = ^/t plus a 
correction of order V~^^^. 



VII. PHASES OF THE LARGE N INFINITE RANGE MODEL 

The infinite range, large A^, model of a spin glass exhibits a transition at a finite (3 from 
paramagnet to spin glass. In the paramagnetic phase, the operator H + \i has a gap. In this 
phase, the behavior is very similar to that of the spherical model considered above. In the 
spin glass phase, H + Xi becomes gapless in the thermodynamic (large V) limit. Assuming k 
components are used to form the ground state of the system, so that the system condenses 
into k states, then the operator H + Xi will have k eigenvalues of order k/V, as well as a 
spectrum of higher eigenvalues. First, will briefly discuss the infinite range model in the 
paramagnetic phase. 

Suppose, as an approximation, that we can look at a system in which Xi is the same for 
all values of i. This is essentially the spherical model considered above. Let this value of Aj 
be denoted A. Then, we find that 

Gu = i{H + Xy% (43) 

Using equation (PSJ), and requiring that Gu = /?, we find the same solution of the self- 



consistency equation as was found in equation (^01). Further, it may be shown that the 
mean-square fluctuations in Gu between different Hij, for this constant A are of order y. 
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Howver, in the paramagnetic phase, it is not true that Aj is a constant. The shght 
fluctuations in Gu away from the average value of equation (p8|) require A, to vary from site 
to site and sample to sample. They produce the difference between the large A^ model and 
the spherical model in this phase. 

We can estimate the fluctuations in Aj. The mean-square fluctuations in Gu from site 
to site and sample to sample are of order 1/V . So, the approximation of taking a constant 
A will almost satisfy the self-consistency equation, but will be slightly in error due to these 
fluctuations in Gu away from the average value. We will flx this by making a small change 
in Aj. Let us consider the change in Gu resulting from a change in \j. This is equal to 

5Gii = Gij5\jGji (44) 

For j y^ i, this vanishes in the thermodynamic limit since GijGji is of order 1/V (this will 
be discussed in the next section). For j = i, this yields 

6Gu = SKP^ (45) 

So, if a constant A produces a small error in the self-consistency equation, we can flx it with 
a small change in Aj on each site. We use equation (|45|) to estimate the change in Aj needed. 
Since the mean-square fluctuations in Gu are of order 1/V, the mean-square fluctuations in 
Aj must also be of order 1/V. The fact that fluctuations in Aj are small, and that fluctuations 
in Gu are small for flxed A, is what permits us to approximate the large N results by the 
spherical results in this phase. 

In the spin glass phase, the spherical model and the large N model differ greatly. In 
the spherical model, the system condenses only into one eigenvalue. The eigenvector for 
this eigenvalue is a randomly chosen vector from a \^- dimensional space of vectors. So, Gu 
can be written as the sum of two parts. The flrst is a contribution from the continuum of 
eigenvalues above the lowest eigenvalue; this part does not fluctuate strongly from sample- 
to-sample (when the spectrum of if + A becomes gapless, this part acquires mean-square 
fluctuations which are much bigger than in the paramagnetic phase, but which are still small 
in the large V limit). The second is a contribution from the lowest eigenvalue. This part 
provides a contribution to Gu of order V{vi)'^, where Vi is the eigenvector corresponding to 
this eigenvalue. This contribution to Gu has fluctuations from sample-to-sample of order 
unity, as (fj)^ has site-to-site and sample-to-sample fluctuations of order 1/V. So, the spin 
glass phase and ground state of the large A^ model will differ greatly from the spherical model. 
The nature of the ground state is the problem that will be addressed with supersymmetry 
techniques later in the paper. It may be assumed that the state in the spin glass phase for 
(3 > '2^/T is very similar to the ground state. 

Using equation (|T^) and equation (|40|), we flnd that the energy of the system in the 
paramagnetic phase is given by 

- vl (46) 



At the phase transition point, this becomes —Vyj^. In the spin glass phase, 
that the energy is given by 
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A ~ a/-, so 




(47) 

The number of spin components A; is in a sense an analogue of Griffiths effects, known 
from finite dimensional systems, in the infinite range model. In the rest of this section we 
will pursue this analogy, to help illustrate the importance of the number k. Consider a finite 
dimensional system of Bose particles with repulsion. As the temperature is lowered, there 
can be a phase transition to a superfiuid phase. If there is disorder, there will also be a 
Griffiths phase, when the system has gapless excitations, but still has no long-range order. 
We can understand the Griffiths phase as follows: some single particle eigenstate of the 
finite dimensional system will be lowest in energy. The particles will begin to condense into 
this eigenstate. However, the eigenstate is localized, and finite in size, so only a microscopic 
number of particles can condense into the state before the interparticle repulsion raises the 
energy of this state, and particles begin to condense into some other state. This leads to 
a large number of states near zero energy, as the states lowest in energy have their energy 
raised by interactions, and produces the gapless spectrum. In the infinite range spin glass 
we consider in this paper, the Griffiths phase is not a separate phase, but the system still 
manages to produce a large number of very low energy states (this number is k) by adjusting 
A.. 

In the spherical model where A is independent of i, the system always condenses macro- 
scopically into the lowest eigenstate of the matrix H, since H + X and H have the same 
eigenvectors. The lowest eigenvalue oi H + X will be of order 1/V, while the next eigenvalue 
will be of order (l/VY^^, as {1/VY^^ is the approximate level spacing in the tail of the 
spectrum of H, so that the system does not macroscopically occupy states other than the 
lowest. 

In the large A^ case, we showed that in the spin glass phase the approximation of using 
constant Aj fails badly. The system system must shift Aj to satisfy the self-consistency 
equation. In the process of shifting, the system will raise the energy of the lowest eigenstate 
oi H + Xi more than it raises the energy of the eigenstates just above it. This means that 
several of the eigenvalues oi H + Xi can become of order 1/V, which leads to A; > 1. This 
process is very similar to the formation of Griffiths states discussed above. 

VIII. CORRELATION FUNCTIONS 

Let us suppose we have a system with a non-degenerate ground state which uses k spin 
components. We would like to investigate the meaning of k and how one can see the effects 
of a given k by looking at the correlation functions of the system. Of course, we have defined 
k as the dimension of the vector space spanned by the vectors 0f . This is some subspace of 
an A^- dimensional vector space. So, if one looks at all V vectors 0i, one can determine what 
k is. Let us instead look at correlation functions of the system. 

First, let us consider the correlation functions within the paramagnetic phase discussed 
in the previous section. Within the large A^ formalism, the correlation function 

T.{m) (48) 
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is equal to 

NGij (49) 

where the average in equation ( ^81 ) is a thermodynamic average for a given Hamiltonian. By 
self-consistency, we have that I](0f0f ) = 1- Now consider Gij for i ^ j, in the paramagnetic 

phase. Clearly, this Green's function vanishes after averaging over different Hamiltonians. 
However, we can compute the mean-square fluctuations in this Green's function, where the 
fluctuations are from sample-to-sample. 

The operator H + Xi has V different eigenvectors. Denote these eigenvectors by ff , where 
/i = 1...V now labels the different eigenvectors. Let each eigenvector have eigenvalue E^. 
Then 

G., = E ^^r< (50) 

Averaging over different matrices Hij will cause the eigenfunctions to have random sign on 
each site, and will cause this quantity to vanish. Further, 

G^,G,. = E ^^^>M (51) 

Averaging over different matrices Hij, we find that the right-hand side of equation ( |5lD will 
vanish, except for the terms in which fi = u. Then we obtain 

G.A-. = E77^(-r<)' (52) 

Since we normalize ff by taking Z](ff)^ = 1, we find that (ff)^ ~ y- So, the right-hand 

i 

side of equation (|5lD contains V terms of order y2, and so GijGji is of order l/V. 

Now, consider the spin glass phase of the system. There will be k eigenvalues which are 
of order k/V. We can again write 

G,G,. = T.T^2i<i^'^^j)" (53) 

as before. However, we can separate out from this sum the terms involving the k eigenvalues 
Efj_ of order k/V. Then, we find that the right-hand side of equation (^) contains k terms 
of order ^. So, GijGji is of order ^. 

IX. SUPERSYMMETRIC FORMALISM 

Here we set up a supersymmetric formalism to determine properties of a large N system, 
averaged over disorder. In the section XI, this formalism will be used to obtain results on 
the value of k for the ground state, the strength of fluctuations in Aj, and other results. In 



this section, we will first derive the formalism in generality for any large A^ system, and then 
discuss simplifications in the low temperature limit. 

To average over quenched disorder, many tricks have been invented in other problems. 
The replica trick and Efetov supersymmetry both attempt to ensure that the partition 
function of a system, for a fixed realization of disorder, is equal to unity. Then, one can 
average over different realizations of the disorder directly. Our goal will be similar. We wish 
to find a method of obtaining A, as a function of H . We will express this using an integral 
over all Aj, such that the integral is equal to unity for any H. 

First we will discuss a naive attempt which does not quite work. Then we will introduce 
a supersymmetry and give the full formalism. 

Consider a given Hij. The goal is to find a set of Aj which satisfies the self-consistency 
equation Gu = (3. One simple way of implementing this might be to integrate over all Aj 
with a set of (5-function constraints on the Green's function. We will also replace equation 
(|n|) with a positivity constraint on G instead. Writing these 5-function constraints as a set 
of integrals over variables Cj, we obtain 

rfA,^e*^»(^--^) (54) 

G>0 

where Ga = (Hy + Aj^jj)^^. However, the above equation is not quite right as there is a 
non-trivial Jacobian coming from the 5-function. To find this Jacobian, assume we have 
found the set of Aj which satisfy the self-consistency equation. Considering a shift 6Xi, we 
would find SGu = Gij6XjGji. 

Therefore, the integral over Aj and c in equation (|^) is equal to l/det(Mjj), where Mjj 
is the matrix given by 

Mij = GijGji (55) 

Note that this is not a matrix square, but a square of individual elements in Gij to produce 
Mij. To cancel this determinant, we can add an additional set of integrals over Grassman 
Oj, Qi variables to equation (B3). This leads to a new equation 



f dXidcidaidaie'^^'^^^^-'^^e^^^^'^^^^^ (56) 

G>0 

The integral in equation ( p6D is equal to unity for any Hamiltonian H. Then, it is 
possible to average over Hamiltonians within the equation. The fact that the integral is 
equal to unity is due to a supersymmetry. There are 2V bosonic variables and 2V fermionic 
variables. Similar supersymmetries have been found useful in other systems [§]. 

It is interesting to compare equation (|56|) to equation (|^). The equations are similar; 
as discussed in reference to equation (^^, in the paramagnetic phase the matrix M can be 
treated as diagonal in the thermodynamic limit. In the spin glass phase the matrix M will 
be very important and will be discussed below and in section XL 

In the infinite range model, we can further simplify equation (|56|). Let us average over 
Hamiltonians H with Gaussian weight e"^'^"'"''^ . We obtain 
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/ dXi dci da, da, dHe-^^^^"^h'''^^^"-f^h''^^^^''^^^' 

J dHe-y-^'^i^' ^ ^ 

For notational simplicity, we will not always write out the factor of r .V^Trg^ • This 

factor will be assumed throughout. At one point in section X, and one point in section XI, 
we will need this factor, but elsewhere it will be left out. 

For given Aj, we can write the integral over H as an integral over O and E with H + X, = 
O^EO (note that the left-hand side is if + A„ not H). Then G = O'^E-^0. Then the 
integral of equation (|57|) can be written 



/ d\i da da, dai dE dOe-^^^'^^'^ J] I^m " ^^1 e^^^^^^'-^^e'^'^'^'^^^^^ (58) 

E>0 ''<'' 

For given k, we will have k of the eigenvalues E^^ scaling as l/jS. In the large f3 limit, 
these eigenvalues are much smaller than all other eigenvalues. It is useful to keep f3 finite, 
but in the large (3 limit we can also assume that Ga is made up only of terms due to the k 
eigenvalues which go to zero, and ignore contributions to Ga from other eigenvalues. Then, 
equation (|^) simplifies. Let us define wf as the eigenvector corresponding to eigenvalue 
E^. In the large /? limit, Ga does not depend on ff for n > k. Further, for /i = l...k and 
p = k + l...f , we have \E^ — E^\ = E^. Also, we can perform the integral over a,, a, to obtain 
det(Mjj), where Mij is defined by equation (|5BD. 

For given k, we will define a quantity Zk- This is what is obtained if equation (^) is 
evaluated in the sector with k eigenvalues tending to zero. Equivalently, Zk is the probability 
that, if we choose H from the ensemble at random, we will have k eigenvalues tending to 
zero. It must be the case that J2k Z^ = I as a. result of the supersymmetry. 

So, in the end we have 

Zk = J dv^dE,S{j:{v^f - l)5(E(t;f<) - l)Z<[v^.E,]Z>[v^,E,] (59) 

i i 

where 

Z^ K, E^] = f dci dai dai H I^m - Ey\ e'^^^^"-^^ (60) 

and 

Z> [t;f , E^] = I dH> dXi e-^-Tr(H2) (det(ii> + x.)ke^^G.,a,G,, (g^^ 

This needs some explanation. The vectors wf are the eigenvectors associated with the 
eigenvalues E'^ for /i = l...k. The eigenvalues E"^ are all of order 1//3. The Green's function 

k 

is given by Ga = J2 ('^j^)^(-^a')~^- The matrix iJ> is a V^-by-V matrix, such that if> + X, 

has only V — k non-zero eigenvalues. The zero eigenvectors of Hy + Aj lie in the subspace 
spanned by the f f . We have H = H^ + J2fi '^{E^Vj. That is, H^ -|- A is what results if one 
looks only at the V — k largest eigenvalues of ii -|- A. 
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Equation (|59|) has been written using two equations: equation ( |60D and equation (^TJ). 
Indeed, equation (|59|) integral almost separates into two parts. For convenience, we will refer 
to Ef^ and ff , for fi = l...k as the lower sector, while we will refer to the part of the integral 
involving if> as the upper sector. Similarly, we may refer to the upper and lower sectors 
of the matrix if + Aj. One way in which the lower and upper sectors communicate is that 
the eigenvectors of the lower sector determine the allowed eigenvectors of the upper sector, 
which, if Aj is not constant, can change the value of Z^. Also, the determinant due to the 
integral over Grassman variables depends on both the upper and lower sector of eigenvalues. 

It is interesting to see how the bound k{k + l)/2 < l^ is realized in this formalism. We 
must show that for k{k + l)/2 > V the integral of equation ( pQ]) vanishes in the large f3 

k 

limit. There are four pieces to this. First, the (5-functions, S{J2 (^f)^(-^/i)^^ ~ P) can be 

k 

rewritten as ^6{ J2 {'^i)'^{PEu)~^ — 1). This then contributes a factor of (4)^, as there are 

V such ^-functions. Remember that E^ scales as 1//?, so fiE^ tends to a constant in the 
large (3 limit. 

Next, the level repulsion between the first k eigenvalues contributes a factor of (■i)'^('^~^)/2. 
The integral over the first k values of £^^ contributes a factor of (4)'^, as Efj_ is of order ^. 
These two combine to produce (i)'=(*^+i)/2. 

Finally, there is the determinant due to the integral over the Grassman variables. This 
is the determinant of the matrix Mij. First consider the part of this matrix due to the 
contribution of the first k eigenvalues to G. This is the determinant of 

E v^v'({E,E^)-'v>;v'; (62) 

This is the sum of k{k + l)/2 distinct matrices, each matrix corresponding to a term in the 
above sum with given ^,v. It can be seen that each of these matrices has one nonvanishing 
eigenvalue which scales as /3^ in the large (3 limit. There is an additional contribution to the 
determinant which is equal to 

2E E v^v'({E,E^)-'v';v^^ (63) 

/i=l ii=k+l 

This is the sum of kiV — k) distinct matrices each with one nonvanishing eigenvalue of 
order (3. There is also a contribution in which both /i and v are greater than fc, but this is 
unimportant. 

The determinant of the matrix Mij is the product of the V eigenvalues of that matrix. 
There are at most k{k + l)/2 eigenvalues which are of order /3^. These arise from equation 

(|6^). The remaining V — k{k + l)/2 eigenvalues are of order /3. So the determinant scales 

as j3y+Kk+i)/2^ 

Combining the various factors of (4)^, (i)'^(*^+^)/2, and ^^+'=('=+i)/2^ -^g ^^id that every- 
thing cancels. However, if k{k + \)/2 > V, there will be a problem. It will not be possible for 
Mij to have k{k + l)/2 eigenvalues of order /3^ since G has at most V eigenvalues. Therefore, 
if k{k + l)/2 > V, equation (|^) will vanish in the large /3 limit. 
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X. A NON-INTERACTING EXAMPLE 

In the previous section, the integral over H + \i was written as an integral over a lower 
sector with eigenvalues of order 1//5, and an upper sector. It will be useful to illustrate 
this technique of dividing an integral over random matrices into two sectors, using a simple 
example. 

Let us consider an integral over all l^-by-V^ matrices if, with weight e"^'^"'"''^ . This 
integral is a Gaussian integral, and can be performed to yield 

We will separate out k eigenvalues from this matrix, breaking the integral into two sectors 
of eigenvalues. Then we will investigate the consequences of requiring that the integral over 



both sectors, including interaction between the sectors, is equal to equation (|6^) . 

Let us pick out k eigenvalues from H, and require them to be equal to E^, /x = l...k. 
Then, the integral over all matrices H can be written as 



ldv^dE>tdH^\{[det{H^ - E,)e-''^^h{Y^{v^f - 1)] l[[\E,- EMT. 



<^J) 



-VtTtHI 



(65) 



where Hy is a (V^ — k)-by-{y — k) matrix. 

If /c = 1, the eigenvalue distribution of if> will be almost unchanged from a Wigner 
semi-circle. Then, the term (det(i7> — E^)Ye~^'^^i^ is approximately equal to 



for all E,, inside the semicircle, — \= < En < -i=- We find 

^^^KY.i^^f - 1) / rfi^>e-^-^^^> (67) 
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/ Vt 2Vt^v-i 



The integral j dvi5{J2i{viY — 1) is approximately equal to 



e^/2(^)V/2 (70) 
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It can be verified that, if tlie last two equations are used in equation ([68|), tliat tlie factors 
of e^ and V^ cancel, as they should. A similar cancellation will be important in the next 
section. 

li k > 1, the term JJ \E^ — E^\ might seem to make equation (|65D greatest when the 

fj,<iy 

eigenvalues i?^ are well separated. However, we know that the original matrix H has V 
eigenvalues, with level density of order V and separation between levels of order 1/V. For 
well separated levels, with \E^ — E^\ >> 1/V, the effect of level repulsion is negligible. The 
solution to this problem lies in the determinants Y\^det{Hy — E^). For a given value of 
E^ the determinant, det(if> — E^j), will slightly alter the level distribution of Hy. This 
will produce a slight decrease in level density around Efj_. This will then, as a result of the 
determinant, det(if> — E^), produce an effective attraction between levels that cancels out 
n \E^ — Ei^\ if \E^ — Ei,\ is much greater than l/V. 

Although this discussion seems elementary, it is worth reviewing this for what will happen 
in the next section. Then, the lower sector of if + Aj will have k eigenvalues all within order 
1/(3 of each other. This will produce a change in the level density of ii>, and will change 
the determinant (det(ii> + Aj))^ away from the naive value one would obtain by assuming 
that Hy has a Wigner semicircle density of eigenvalues. The true level density of Hy will 
be evaluated in the next section and in the appendices. Since k will be much greater than 
one, this change in level density will be very important. 



XI. GROUND STATE FOR THE LARGE N INFINITE RANGE MODEL 

Having set up the supersymmetric formalism in the previous section, we will use this 
formalism to determine the most likely value of k and the strength of fluctuations in Aj. It 
will be shown that the system will obey k oc V'^^^. Further, it will be shown that the mean 
square fluctuations in Aj are equal to -^. We will also derive a slight change in the average 

value of A from J^, and a slight gap in the low energy excitations. 

The starting point for the calculation is equation (p^l). We will evaluate this integral 
for different k and use this to determine the most likely k. For any given k, we can write 
the value of this integral as Z^ = e^'=. Below we will evaluate various contributions to 
Sk- Due to the complexity of the calculations, we will only evaluate terms in Sk to order 
kV, V^logV, Vlogk, V, k'^logV, and k'^logk. There will also be terms in Sk of order k"^, V/y/k, 
and smaller, which we will not evaluate. All the terms of order kV, yiogk, V^logV, and V will 
cancel. It will turn out that the other terms do not cancel; given these terms, we will show 
at the end of the section that k is of order V'^^^, with the exact ratio k/V"^/^ undetermined 
without a computation of terms of order k"^ and V/yk and in Sk- Note that for k of order 
V"^/^, then k"^ and V/yk are of the same order. 

The procedure to evaluate equation (|59D will be to start with all Aj equal to each other. 
We will then write Aj = A. In this case, we will evaluate the integral over the various 
eigenvalues E and eigenvectors v to obtain the result. Then we will consider the effect of 
site-to-site fluctuations in Aj at the end. From this we will get an effective action for A which 
will tell us the average value of Aj and the strength of fluctuations about this average. 

The rest of this section will be divided into four subsections. First there will be a 
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calculation of the contribution to Zk from the upper sector. Then there will be a calculation 
of the contribution to Z^ from the lower sector. These two calculations will be for the case 
with all Aj equal to a constant A. In the third subsection, we will consider the effect of 
permitting A, to vary from site to site. In the final subsection, we will put the calculations 
of the first three subsections together and obtain final results. 

A. Upper Sector 

For constant Aj = A, Z^ is completely independent of the eigenvectors vf for /i = k+l...V . 
Zk does still depend on E^ for fi = k + 1...V. The matrix H + X has k eigenvectors near zero 
and V — k other eigenvectors which are all positive. The average value of A is determined 
by a balance of two effects. First, there is a term in the probability distribution for A (when 
all Aj are set equal to A) like e~''^'^^ . This term simply arises from the requirement that H 
must have k eigenvectors equal to —A and from the term in equation (pT]), e~^'^'^^^ , where E 
are the eigenvalues of H. This term favors a small value of A. There is also a term due to the 
integral over E^^ for fi = k + 1...V. The requirement that these eigenvalues all be positive, 
and the effect of level repulsion between these eigenvalues and the first k eigenvalues, tend 
to favor a larger value of A. 

For example, in the simple non-interacting problem of the previous section, when we 
removed one eigenvalue, the level repulsion and the Gaussian weight canceled, and the 
eigenvalue could have been found anywhere between —J^ and J^. In the case considered 
in this section, the effects of level repulsion will be stronger, as k eigenvalues will coalesce. 
Further, the positivity constraint must be satisfied. We will also find from this a slight 
difference between A for the large N model and A for the spherical model, which will imply 
a shift in energy from equation (|12D- 

The integral over the upper sector of eigenvalues is performed in appendix A. There, 
it is shown that the eigenvalue distribution of the upper sector is slightly distorted from 
the Wigner semicircle, and that A is slightly shifted from ^ to ^ — |^2^/^(p:)^/^. If the 
eigenvalue distribution of the upper sector were unchanged from a Wigner semicircle and 
A = J^, the determinant (det(if> + A))'^ would equal 

/ ^ \k(V-k)^k(V-k)/2 /y-^N 



Also, e~^^'^^ would equal e~^^ . However, the slight change in the eigenvalue spectrum, and 
slight shift in A lead to 



(det(i7> + \)fe-'''^'' = ^^)Hv-k)^kiv-k)/2^-kv^k^-l^ ^^2) 



as can be shown using the results in appendix A and calculating the determinant by inte- 
grating the eigenvalue density from appendix A. 

In appendix A we calculate A and a, where a is the lowest eigenvalue of /f>. The sum, 
A-|-a gives a gap to excitations. The gap is equal to ^2^/^-4^ (■^)^/^ and will vanish as a power 
law in V . In figure (||) we illustrate the changed eigenvalue distribution, superimposing the 
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Wigner semicircle of the original matrix H for comparison. The gaps in the figure are 
exaggerated and are not to scale; the figure is approximate only. 

The factor of | in the various above equations should not be a surprise. For a matrix 
with a Wigner semicircle distribution of eigenvalues, there is a square-root singularity in the 
eigenvalue density near the end of the semicircle. This means that the lowest k eigenvalues 
range in energy up to a distance of order [k/VY^^ from the end of the semicircle. Having 
changed the distribution of eigenvalues from the Wigner semicircle so that the lowest k 
eigenvalues all lie at zero energy, the eigenvalue density of the rest of the matrix will adjust, 
but once one looks at energies several times {k/VY/'^ from the tail, the change in the 
eigenvalue distribution is small as a result of screening. This argument gives the scaling 
of the gap. 

To understand the difference between equation (|7ID and equation (|72[), we must calculate 
g-feVrA (^({q\^(^]^^ _|_ \Y^ for given A. The naive evaluation of this is what gave equation (^TJ). 
The evaluation of the determinant is very much like evaluating the potential energy of a 
system of charges interacting with a logarithmic Coulomb interaction. We can then write the 
true eigenvalue distribution of if> as the sum of two pieces: a Wigner semicircle extending 
from — Y ^ to w^, and a negative term extending to a distance of order [k/VY^^^ with the 
integral of the negative term being k. The negative term can be thought of as a "hole" . 

The interaction with the Wigner semicircle, multiplied by e^^'^^ , is independent of A, 
for A inside the semicircle, and the product is equal to equation (ffll). So, we can think of 
equation ( [72|) as arising from equation ([7l|) , multiplied by a contribution arising from the 
interaction of the k levels of H at —A with the "hole". This gives equation ([721). 

Finally, since we are considering the integral over iJ>, recall that in equation (^9t), we 
have normalized the integral by 

(73) 




In the calculation of equation ( [7^ ) and of appendix A, when we integrate over Hy, we gain 
an extra factor of 

|rfi7>e-^^^(^') (74) 

The ratio of these produces a factor of 

(75) 

TT 

This is similar to the calculation in section X. 

B. Lower Sector 

The rest of equation (|59D involves the lower sector of eigenvalues. Let us write 

(76) 
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Here, 0f is a ground state of the system, up to a change in normahzation (compare to 
equation ([T5|) . Further, let us change variables and write i?^ = p~^. Then, the lower sector 
of the integral can be written as 



^1<U 



(77) 



If we ignore the determinant in equation ([77|), the integral over 0, at fixed p, can be 
performed using a series of Lagrange multipliers to implement the 6 functions. We can take 
the integral 

' 0rn^(E(</>r)'-/?)^(E(</>n'-p.))n(^(E(</'rc)) (78) 



and write it as 



!■ 



271 271 271 



We can integrate over the fields 0f , to obtain an action for the Lagrange multiplier fields 
c. The integral over c can then be handled by the saddle point method. One possible saddle 
point is Cfiu = 0, c^ = iE^V/2, Ci = 0. There is a zero mode for fluctuations about this saddle 
point which enforces the constraint J2i_i Pn = 1^^ ■ This zero mode is a mode in which all the 
Ci are increased by some constant, while all the c^ are decreased by the same constant so 
that c^ + Ci is unchanged for all p and i. 

At the saddle point, we find that the integral given by equation (|79D is equal to 



e 



"'=n(T)' (80) 



For V much greater than fc, we can look at other fluctuations about the saddle point in 
a Gaussian approximation. The fluctuations in c^p contribute a factor of 



r dCf^ ^-^c%p^p^ ^ ^'^yk{k-i)/4(^2n)-'''^k-i)/2 [](p^/\/)-('=-i)/2 (81) 



The fluctuations in c,- contribute a factor of 



The fluctuations in c^ can be ignored, as there are only k such terms fluctuating, and this 
would only lead to corrections to 5*^ of order k. We are not considering terms in 5*^ that 
small. 

We will now evaluate the contribution to Z^ from det(Mjj). This determinant involves 
both the lower and upper sector of eigenvalues. As discussed in section IX, this matrix has 
k{k + l)/2 eigenvalues of order /3^ and {V — k{k + l)/2) eigenvalues of order [3. The trace of 

26 



y/32 



the matrix is VjS'^. So, the k{k + l)/2 largest eigenvalues can be at most of order 

Let us assume that in fact each eigenvalue is approximately -p-; making the approximation 
only leads to errors of order k"^ in S^- 

Since we will find in the end that k oc V"^^^, we have V — k{k + l)/2 >> k{k + l)/2. So, 
to find the product of the other V — k{k + l)/2 eigenvalues, we can simply throw out terms 
in Mij of order /?^ (these are terms shown in equation (|62l)) and calculate the determinant 
of the remaining V-hy-V matrix. We find 

G.,G,. = 2X: E -^<^K^-=2E E ^'Pimv] (83) 

^=1 a=k+l -^M-^a 11=1 a=k+l ^<^ 

Here, E^ is of order 1//5, while Ea is from the upper sector of eigenvalues. Averaging over 
0f , we can ignore terms in which i ^ j and obtain 

GuGu^2 E /5-^<^J (84) 

a=k+l ^o- 

This is a Green's function of the matrix H^. The Green's Gu function oi Hy is approximately 
2i/r and so equation (|^) is equal to 

4/3v^ (85) 

Putting these together, we obtain 

det(M,,) ^ { ^Jf^.^^/ ^'^''>/\(34V^r-'^'^'^/' (86) 

There will be corrections to this which lead to corrections to Sk of order V/vk. These 
corrections will be discussed at the end of this section and in appendix B. 

We now need to know what the spectrum of E^ is, for /x = l...k, to determine what 
some of the contributions calculated above are equal to. The exact spectrum is difficult 
to determine, but we can determine enough for our purposes. There is an effect of level 
repulsion which tends to push the energies apart, but it may be seen that the other terms 
in equation ( [77|) are greatest when the energies are all equal to k/{VP). The effects of level 
repulsion will push the E^ apart, but on the other hand if we considered only the effects 
of level repulsion, ignoring all other effects except the requirement that all eigenvalues be 
positive, we still want all E^j, to be of order k/(yp), up to factors of order unity. So, 
considering all effects, we must have that £^^ is of order k/(yp). 

Also, we would like to know how large \E^ — E^\ is. Although each E^ is of order k/{VP), 
the separation between energy levels may be much smaller, as the contribution of equation 
(^) depends on E^ and is greatest when all E^ are equal. There are other terms which 
have a similar effect and are also greatest when all E^ are equal. However, we will now show 
that despite this |i?^ — Ei,\ is still of order k/{VP). It can be shown that, if equation ( |52D 
is expanded out in SE^ = E^ — k/{Vl3), that there will be term of the form 



4^ yav ! 



w> (87) 
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Thus, the energy levels E^ repel each other due to the absolute value, \E^ — Ey\, but are 
confined by a potential that, for E^ approximately equal to k/(y(3), is given by equation 
(p7|). If we have k levels, E^, repelling each other, in a confining potential given by equation 
(^), then for k of order V^^^, or greater, we find that the separation between levels is of 
order k/Vp. So, we can assume that \E^ — Et,\ is of order k/{VP) since we will find in the 
end that k = l/^/s > yi/s^ 

Putting all this together, we obtain a contribution to equation ( |5UD from the lower sector 
given by equation (|77D and a contribution from the upper sector given by equations ( [72| , [f5D , 



as well as the determinant of equation (^). All this is at fixed A. Equation ([7^) includes 
a number of terms due to the change of variables and from the level repulsion, as well as 
contributions evaluated in equations (|3,^,|^). If we take typical p^ to be of order V [3/k 
as discussed in the above paragraphs on level repulsion, we find 

Zu^{'^f'\^r'^'^{^Y'' (88) 

k'^ V IT 



C. Site to Site Fluctuations in Aj 

The above calculation was performed for given A. Finally, we must obtain the strength 
of fiuctuations in A, and perform the integration over A. We will now obtain an effective 
action for A. Let Aj = A + SXi, where A is the average value of Ai calculated above. Then, 
J2i <^Aj = 0; this is important and will be discussed more below. The action for Aj arises 
from the change in (det(i7> + A))'^, as well as the change in H required to ensure that H + X 
still has k zero eigenvalues. 

First, let us evaluate the change resulting from the determinant. We have 

(det(i7> + A))'= = e^T^i°s(^>+^») (89) 

Expanding the log in the above equation in 6Xi we obtain to second order 

gfcTrlog(//>+A)g-fciG2^5A2 ^ ^feTr log(//>+A)g--2fcr5A2 /ggN 

Finally, remember that H + Xi must have k zero eigenvalues. For Aj = A, a constant, 
this contributed a term e~'^^'^^ as found in equation ([72|) . Let us write H in a. basis of 
the eigenvectors of if + Aj. Then, since {H + Aj)f'^ = 0, for /i = l...k, we must have 
Hfi^a = Hi'^i^i^i, for /i = l...k and a = 1...V. Also Ha,^ = J^iViXiV^. Then, we must 
evaluate 

-Vr f: E (Hl^+Hl,) -2Vr f] E £. -^^^<? 

g /i — la — 1 —- p /^— 1 a — 1 fQl") 

Since (ff)^ and (f")^ are approximately 1/V , this is 

Multiplying equations (0,^), and integrating over Aj, we get find an extra contribution 
to equation (p^) 
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'Akr' 



(93) 



Corrections to this will lead to corrections to 5*^ of order V jyk. These will be discussed 
later. 

The result from equations ( p0ip2|) is interesting in that it requires that Yl,i^\ = 0. If 
we evaluated the shift in determinant for constant shift in A (all 6Xi equal to the same 
constant), it would exactly cancel the change in e~^'^^ . However, here we would find that 
the perturbative method used to obtain equation (0) was not valid, and that constant shift 



in A could violate equation (ITJ)- Further, we would not have the factor of two in equation (92) 
for constant shift in A. The factor of two arises since the shift due to 6Xi involves off-diagonal 
terms in H, while a constant shift in A would involve diagonal terms in H. Effectively, if 
we look at fluctuations in Aj, the mode in which all Aj fluctuate together is very different 
from the other modes in which J2i ^^i = 0. The mode in which all Aj fluctuate together will 
be discussed in the next section; this mode is important for considering fluctuations in total 
energy. 

D. Results 

The final result, combining equations (^) and (PBj), is that Zk is of order 

In obtaining this result a large number of terms have canceled. It may be verified that all 
terms with an exponent of order kV or V cancel. Further, it may be verified that r cancels 
exactly, as required. 

Equation (0) is maximized for 

k oc y2/5 (95) 

There will also be terms contributing to Zk which look like e*^^ , for some constant a, as well 
as exponential of lower powers of k'^ and powers of V. To evaluate these terms is beyond 
the scope of the present calculation. These terms can shift the constant of proportionality 
in equation (pSj), but cannot change the power law dependence of k on V. It would be nice 
if a more detailed calculation could obtain the exact proportionality constant, as well as 
checking that J2k ^fc = 1- The requirement that J2k ^fc = 1 can be verified here only so far 
as noting that for k oc V'^^^ that Z^ does not have any terms in it like e'^ '"s'^ or like e^, but 
we cannot check for cancellation of contributions to 5*^ of order fc^ = V^^^ or lower order. 

Further, there will be contributions to Sk of order V/vk that would arise from a more 
careful calculation of det(Mjj) and from a better calculation of the integral over fluctuations 
in Aj. We will discuss these terms in appendix B, and show that they are in fact of order 
V jyk. For k = V"^^^, these terms will be of the same order as the /c^ contributions to Sk- 
For k < V"^^^, these terms will be the most important contribution to 5*^, and since they 
are negative, we will find that 5*^ is negative and Zk is much less than 1; for k > f^/^, 
equation (Q) will give the most important contribution and will also lead to Zk being 
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small. Therefore, although terms of order V/yk can change the proportionality constant 
in equation (0), they cannot change the power law dependence of k on V . These terms 
can also lead to a change in the gap, A + a, derived in appendix A, but will not change the 
scahng of the gap with V . 

The results on the ground state can be summarized by equation (^) for the value of /c, 
equation (^) for the fluctuations in Aj, and the calculation in appendix A for the average 
value of A, and the weak gap to excitations above the ground state. 

XII. DISCUSSION OF RESULTS ON THE GROUND STATE AND 
CONNECTION TO PROPERTIES IN THE SPIN GLASS PHASE 

The two most important results are the value of k and the strength of fluctuations in 
A. We will discuss these results in this section and try to interpret them. Next we will 
discuss fluctuations in the total energy of the system, as discussed at the end of section XI, 
subsection C The shift in A is also interesting all will be discussed next. Finally, we will 
argue that the results of the previous section are also applicable to the spin glass phase as 
well as the ground state. 

The value of k is perhaps somewhat suprising. One might have expect k oc yV as we 
have the bound that k{k + l)/2 < V . However, imagine considering a different problem from 
that of finding ground states of the system. Suppose we had instead looked for stationary 
states. These are configurations of the spins such that the energy does not change for small 
changes in spin configuration. The problem of finding stationary states is the same as the 
problem of finding ground states, except that equation ( ]TT| ) does not need to be satisfied 
any more. Then, consider a solution in which all Aj are equal to some constant A. Now, we 
can take A to lie in the middle of the Wigner semicircle, instead of near the edge. Then, the 
solution of the problem would proceed exactly as before, except that the corrections to the 
determinant calculated in appendix A would be of order 

{\f'' (96) 

without the factor of 2/3 in the exponent (compare to the discussion of the origin of the 
factor 2/3, as this factor arises near the edge of the semicircle). In this case, it may be 
verified that Zk is greatest for fc^ of order V . 

Since the same bound on k holds for stationary states as for ground states, the most 
likely value of k for stationary states is not surprising, as here k"^ is of order V . However, 
one can perhaps argue that typical k should be less for a ground state than for a stationary 
state. From most stationary states one can reduce the energy by aligning some spins and 
reducing fc, bringing one closer to the ground state. This argument is of course very crude, 
but it is an attempt to interpret the result for k, especially since the calculation of section XI 
and appendix A shows that the result for k depends crucially on the square-root singularity 
of the density of states of H. 

In fact, it is not surprising that the square-root singularity should lead to a reduction in 
k. We expect that A will be near J^. The system must take k eigenvalues near —J\ and 
bring them to zero energy, by tuning Aj. The square-root singularity implies that there is 
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a lower density of eigenvalues near the end of the semi-circle, and so it is more difficult to 
bring k eigenvalues to zero energy. When looking for stationary states, instead of ground 
states, we can take A near the middle of the semicircle, and there are many more eigenvalues 
around to bring to zero energy, so k can be bigger. 

More interesting is the strength of fluctuations in A. We have shown that the mean square 
fluctuations are of order 1/k. The fluctuations in A on a given site measure the fluctuations 
in the energy of the bonds connecting that site to other sites. Each site is connected to 
V other sites. The term in the Hamiltonian connecting the two sites has mean square of 
order 1/V. However, the energy due to a bond connecting two sites is determined both by 
the strength of the bond and the correlation between the two spins on the sites. Above, 
we estimated that the mean square correlation function between two sites is of order 1/k. 
Putting all this together, the energy on a given site is the sum of V terms, each with mean 
square fluctuations of order y^. So, it is not surprising that the mean square fluctuation in 
A for a given site would be of order 1/k. 

The total ground state energy of the system is given by the sum of the different Aj. Since 
J2i^K = 0, we must look at fluctuations in A to get the fluctuations in the ground state 
energy. From the calculation of appendix A, and from the discussion of the "hole" in section 
XI, we understand that A is being attracted to this "hole" which has width (y)"^'^ = k~^. If 

A increases beyond J^ then e~''^'^^ (det(if> + A))'^ becomes small since the contribution to 
the determinant from the semicircle decreases rapidly for A > J^. On the other hand, if A 
gets much smaller than the value calculated in appendix A, A = — -\/t + |a/t2^/^/c^/'^, it will 
be difficult to satisfy equation (11). So, A varies over an interval of range k~^. As A is varied 



over this interval, one must find how e~^^'^^ (det(if> + A))*^ varies; this can be accomplished 
by a calculation along the lines of that in appendix A, or by using the following simple 
argument. The attraction between A and the "hole" is effectively a logarithmic Coulomb 
attraction between two charges of opposite sign and magnitude k. The hole has width k~^. 
A simple calculation on this Coulomb system shows that A can vary only over a region of 
width k'"^, where one factor of k~^ arises from the width of the "hole" and the other factor 
of k~^ arises from the magnitude of the charges. So, the fluctuation in the energy of the 
system is proportional to Vk~^ = V^'^ and the square fluctuation in energy is proportional 
to V'/\ 

This may at flrst seem surprising. For flnite dimensional random systems, the square 
fluctuation in the energy is usually of order V by the central limit theorem; the total energy 
is the sum of V different quantities with fluctuations of order unity. However, in an infinite 
range model such as is considered here, each site is connected to V different sites, and the 
fluctuations in the energy for a given site are correspondingly much smaller. Therefore, it 
is possible for the mean square fluctuation in total energy to be of order V'^^^ instead of V. 

In the paramagnetic phase, where the mean-square fluctuations in A^ are of order 1/V, 
the mean-square fluctuation in total energy is of order unity. 

The shift in A, leading to A slightly greater than —J^ is not surprising. This implies that 
the ground state energy of the large N system is slightly higher than that of the spherical 
model. This is of course not surprising since the spherical model has additional freedom to 
choose its ground state, as the constraint on the sum of the spins from the spherical model 
is less constricting than the large A^ constraint on the length of each individual spin. 
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Having obtained these results for the zero temperature problem, we expect that they will 
hold true for the spin glass phase. As discussed above, the thermodynamic limit and large 
(5 limit both involve sending eigenvalues to zero, so the properties of the spin glass phase 
should be very similar to the zero temperature problem. 

XIII. CONNECTION WITH EFETOV SUPERSYMMETRY 

We will discuss a slight modification to the previously constructed supersymmetric for- 
malism. This modification will make the formalism very closely related to the Efetov super- 
symmetry used in calculations on non-interacting disordered systems. It opens the possibility 
of using field theory techniques to deal with disordered and glassy A^ systems in the large N 
limit; this might prove much more useful for finite dimensional systems when the techniques 
of random matrix theory are not available. Unfortunately, in the spin glass phase of the 
infinite range model it has not yet been possible to proceed with this modified formalism 
to the extent that has been done above with the original formalism. Still, the technique is 
interesting in itself, and will provide a field-theoretic method of obtaining the properties of 
the paramagnetic phase discussed above. 

One complication when trying to apply field theory techniques to equation (^) is the 
appearance of the Green's function in the exponential. This makes it impossible to do what 
one would normally like to do, namely integrate over H before integrating over any other 
variables, as it is difficult to integrate over matrices Hij when the action depends on the 
inverse of the matrix. We can rewrite equation (|56D in a slightly simpler form that avoids 
some of these complications. First we can rewrite it as 



j dGijd\ — dUi da,e''^^^^^-'^e''^^^^'''^^^6{G - {H + X)-') (97) 

J 2ti 



We can replace the (5-function by S{{H + X)G — l)det{H)^ . Writing the (5-function by an 
integral over a set of auxiliary variables Cij and writing the determinant as an integral over 
a set of Grassman variables ipij,ipij, we obtain 

/ dG,jdCijdi^ijdij,jdXi dCi dUi t/a.e*'='(G,,-l)ga,G.,.a,G,,g*TV{c(HG-l))gTr(^H^) ^gg^ 

G>0 

This introduces a second supersymmetry into the problem, connecting the variables c, G 
with the variables ip, ip. This supersymmetry is closely related to the Efetov supersymmetry 
used in calculations for disordered system [Q. 

In the paramagnetic phase, it is possible to use this formalism very simply. It has not 
yet been possible to obtain anything interesting in the spin-glass phase using this formalism. 
In the paramagnetic phase, we can integrate over if, and then introduce a supermatrix field 
Q to decouple the integral, following procedures similar to those used in the Efetov super- 
symmetry technique 0. In the paramagnetic phase, a simple saddle point approximation 
on Q suffices, since Q has a gap for fluctuations. At this point, the notation becomes quite 
complicated, although the ideas are simple, and not very different from those used in the 
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Efetov technique used for non-interacting systems. We will simply state that such a saddle- 
point technique reproduces all the results of the section on the paramagnetic phase of the 
large A^ spin glass. 

In the spin-glass phase, the fluctuations of the superfield Q become gapless, and it 
becomes difficult to proceed with this technique. We expect that the bosonic sector of Q 
may have k eigenvalues separate out from the others; there will be then V — k remaining 
bosonic components and V remaining fermionic ones. The difference of k between these may 
produce the factor of (det(if> -|- Aj))^ that was so crucial before. However, it is difficult to 
use Efetov supersymmetry near the tail of the spectrum of eigenvalues of a random matrix. 
This is the problem in proceeding further in the formalism of this section. 

XIV. TROUBLES WITH REPLICA SOLUTIONS TO THE PROBLEM 

We will briefly review the replica solution to the large N spin glass problem [^, indi- 
cating a number of mathematical problems. The replica solution yields good results in the 
paramagnetic phase, but does not correctly address the spin glass phase. Two important 
questions considered in this paper, namely the number of spin components used to form the 
spin glass state and the site-to-site or sample-to-sample fluctuations in A, are not even con- 
sidered within this replica formalism. However, we can point out a few other mathematical 
problems. 

In the replica technique we compute a partition function such as 

Z= I #f'"e^^ (99) 

where (j)f has now been given an extra index a. This is a replica index ranging from l...n. 
In the end, the limit n = is taken. 

After averaging over H, one can then decouple the average with a replica matrix Qa,f3,ii,u- 
The result is 

Z= I t/0r^Q«,/3e-''^^^(^')e<^"'^«-''-'"^"' (100) 

Finally, the constraint on the length of the spins is enforced with a Lagrange multiplier A"" 
to obtain 

Z= J dX'^d<prdQa,pe-''^^<'^"h'f'r(Q^^^-.-^+^^m\-^^>^? (101) 

Taking saddle points in both Q and A, one can obtain results in the paramagnetic phase. 
However, there are several problems. First of all, the decoupling matrix has a number of 
components of order A^^. The fact that Q has a gap to fluctuations in the paramagnetic 
phase of order V is not very useful in the large N limit, since the number of components is 
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so great. It is interesting to note that restricting A^ to be of order \/V ^ remembering the 
result on the maximum k to form a ground state, we would find that the matrix Q must 
still have of order V different components. 

The fluctuations in Q do have at least one important effect. If one ignored these fluctua- 
tions, one would naively think that the gap for fluctuations in Aj is of order A^, since for flxed 
Q this is indeed the correct result. However, integrating over the A^^ different components 
of Q must reduce the gap for fluctuations in Aj to be of order V ^ as this is the result found 
in section VII. 

Another problem has to do with the spin glass phase. In the supersymmetric formalism, 
we had to restrict to the sector H + \i > 0. This might seem to be a problem with the 
supersymmetric formalism, as this is a slightly strange restriction to enforce. However, a 
similar problem will arise in the replica formalism. The procedure used is to integrate over 
all H and then decouple this to produce an integral over Q. Then, an effective action for 
Q and A is obtained by integrating over 0. However, since the flelds are bosonic, the 
integral over is only well deflned if iJ + A > 0. So, the procedure of integrating over H, 
decoupling, and then integrating over 0, is ill-deflned, unless one also can somehow apply 
the same restriction if + Aj > 0. One cannot integrate over all H at flxed A and still be 
assured that H + X will be positive deflnite. 

Further, in the spin glass phase, we saw from the previous sections that H + X is close to 
gap less; the gap vanishes as a power law in the thermodynamic limit. The replica matrices 
also will have gapless excitations. Considering the calculation of appendix A, we expect 
that a detailed consideration of fluctuations is crucial to getting correct answers here, and 
cannot easily be obtained by using replica techniques, or techniques outlined in the section 
XIII. 

Another problem lies in the order of limits. There are the limits n ^ 0, A^ ^ oo, and 
y ^ oo. Even without the replica limit, there are complications in interchanging the large 
A^ and large V limits, as we have found in this paper. 

In conclusion, although the replica technique is very powerful for certain problems, es- 
pecially problems with interaction and flnite A^, where other techniques cannot be used, the 
old results in the literature for the large N spin glass using replica techniques should not be 
trusted fully. Although very useful and enlightening in the paramagnetic phase, in the spin 
glass phase there will be problems. 

XV. HERMITIAN AND SYMPLECTIC SYSTEMS 

Instead of the system of real Hamiltonians considered in above, we can consider more 
general systems. One possibility it to take the spins to be A^ component complex vectors, 
instead of A^ component real vectors. Then, we can take the matrix Hij to be an arbitrary 
Hermitian matrix instead of simply a real symmetric matrix. So, the system will be deflned 
by the Hamiltonian 

H = J2^':<P^H,, (102) 

where i,j index the various sites and range from 1 to V, while /i indexes various components 
of the spin and ranges from 1 to A^. We have the constraint that 
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Y.m' = N (103) 

We will simply sketch the changes in the Hermitian case from the real case considered 
above. Previously, we derived the bound that k{k + l)/2 < V where k was the number of 
spin components needed to form a ground state for a given H. This was based on the fact 
that for a real matrix iJ + A, to obtain k zero eigenvalues requires k{k+l)/2 free parameters. 
For a general Hermitian matrix H + X, obtaining k zero eigenvalues requires k"^ parameters. 
So, we may obtain the bound that k'^ < V. Here, k counts the number of complex spin 
components. 

Similarly, a supersymmetric formalism may be derived for the Hermitian case. The 
formalism of section IX requires works equally well for the Hermitian case. The only change 
required is that we will write H + X = U'^EU, instead oi H + X = O^EO. Then, the 
measure of integration will be different in the Hermitian case, so that we will write H \Efj.— 

E^P|det(//> + A)p'=, instead of n \E^, - ^^||det(iJ> + A)|'=. The calculation of section XI 

can still be pursued for the Hermitian case, and the power law dependences will be the same 
in the Hermitian case as in the real case considered above. 

The formalism of section XIII must be slightly changed for the Hermitian case. The 
matrices G and H can now be complex. To implement the constraint 6{HG — 1) we will 
then need to use a complex field c. Using complex G and c amounts to doubling the number 
of bosonic degrees of freedom. We will then need to use two fields ip and two fields ip to get 
the correct Jacobian. However, all these changes are easy to implement. 

There is also a symplectic case, in which H is taken from the symplectic ensemble of 
matrices. In the symplectic case, we have the bound 2/c^ — k < V. Also, the level repulsion 
factor in section IX must be replaced with H \Eij, — E^\^\det{Hy + X)\'^'^. Similarly, compared 

to the real case, the symplectic case will require four times as many fields using the formalism 
of section XIII. These changes are also all easy to implement. 



XVI. CONCLUSION 

In conclusion, we have looked at the problem of the large N infinite range spin glass, and 
obtained several results on the spin glass phase. One reason for the interest in the problem 
is that it is perhaps the simplest system which combines both disorder and interaction. It 
is simple enough that one feels that there should exist a simple solution, but, as discussed 
above, previous attempts to solve this problem do not fully succeed. It is perhaps surprising 
that such a simple problem requires a technique as complicated as that used here. 

This problem exhibits many interesting aspects of disordered systems. It has an analogue 
of Griffiths effects and a non-trivial ground state. The infinite range spin glass with finite 
A^ has been solved by replica techniques before. However, it is nice to have a glassy sytem 
that can be solved without using replica techniques, to use as a check. For certain aspects of 
random matrix theory, such as level-level correlation functions, replica techniques run into 
troubles, so it is interesting to find a supersymmetric technique to deal with a glassy system. 

The results derived include the scaling of k as V"^^^; the scaling of the gap to excitations 
as (^)^^^ = V~'^^^; and the magnitude of fluctuations in Aj. It has not yet been possible to 
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derive the exact ratio of k/V'^^^ . It is interesting that we deal with many quantities which 
vanish in the thermodynamic hmit in this paper. The gap, the ratio k/V , the correlation 
functions in the spin glass phase, and other quantities all vanish in the large V limit, but 
they do so more slowly than 1/V , which makes them interesting. 

The method introduced in the paper will hopefully be of more general use. It would be 
very interesting to try to use this technique to attack problems with finite range interactions. 
Finite dimensional, large N spin glass systems are one possibility. Another possibility are 
finite dimensional, large N system with disorder, but without frustration. One example of 
these system, the large N dirty boson model Q], has been treated with an RG technique. 
The supersymmetric technique outlined here might be useful in that context. It would also 
be very interesting if a method were found for including 1/A^ corrections in this formalism. 
For finite N in the large V limit we expect replica symmetry breaking, so it is not clear 
what would happen given the first 1/A^ correction to the large N limit. The self-consistency 
equation could have multiple solutions for finite iV, and it is not clear what would happen. 

XVII. APPENDIX A: SOLUTION OF VARIATIONAL EQUATION FOR A 

Here we consider a problem in which all Aj are equal to each other. We will write Aj = A. 
The problem is to find the most likely value of A. The reason for considering this problem 
is discussed at the start of section XL 

Equation (^) depends on A in two ways. First, there is the factor of e~^'^"^^ . Second, 
there is a factor of (det(iJ> + A))'^. As discussed at the start of section XI, we must compute 

e-^^"^'(det(i7> + A))^ (104) 

Let if> have a density of eigenvalues equal to V p{x). Since if> has V — k eigenvalues, 
we will have J p{x)dx = ^^y^. The problem of finding the density of eigenvalues of if> and 
the correct value of A reduces to maximizing the functional 

-VkrX'^ + V^- dx dy p{x)p{y)\og\x - y\ + V'^ dx p{x){-Tx'^ + — log|x + A|) (105) 

where the first integral represent the effects of level repulsion and the second integral repre- 
sents the effects of the Gaussian confining potential and the determinant (det(iJ> + A))'^. 

First, we will find the density of eigenvalues of if> for given A. Then, we will find the 
optimal value of A. By varying equation ( |105| ) with respect to p(x), and then differentiating 
the resulting functional equation with respect to x, we obtain 

J dyp{y)log\x-y\=2TX-——— (106) 

Equation ( |106| ) only holds for x such that p{x) > 0. We assume that p{x) is non-zero for 
a < X < b, and zero for a; < a or a; > 6. 



The techniques for solving this equation can be found in the book ||T0| . We simply quote 
and use the results of that work here. 

The function p{x) can be obtained from the expression 
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(107) 

This integral can be performed by a series of simple steps: write y = ^-^z + ^^ so that the 
integration range for z extends from —1 to 1. Then write z = cos{6) to get an integral over 
6 from to TT. Finally write w = e'^, for a complex variable w and extend the integration 
range for 6 from 0...tt to 0...27T. Then, the integral over w is an integral over the unit circle 
in the complex plane and can be performed by contour integration. 
The result is 

p^jr) = i^(x-a)(6-a;)(2r + , \^^^ .] ) (108) 

For k = 0, this reduces to the well known Wigner semi-circle. 

There is also a consistency equation that must be satisfied that helps determine the 
limits a and b. This equation is 

(109) 

This integral can be performed by the same contour integration techniques as the previous 
one. The result is 

= r(a + b) , (110) 

VJ{a + X){b + X) 
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In the case A; = 0, this reduces to the requirement that a + 6 = 0. 

b 

nt nr\ nnr — - 

a 

using contour integration. The result here is 



Finally, we have the constraint that / p{x) dx = ^y^- This integral can also be performed 



^-^ ri'-^r-k + ^^±^ (111) 

\J{b + X){a + X) 



V ^ 2 ' V 



For A; = 0, this equation gives b — a = 2~7^\J^^^- 

Putting equations ( [110|jlll| ) together, we can obtain a and b. Let us write 



2 + / , , 

b-a = ^^ (112) 

X = -a + ^ (113) 

We assume that / and 5 are small numbers. For A; = 0, we will find that / is of order k/V. 
We approximate J{b + A)(a + A) = -^. Then equation (|110|) becomes 
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Ignoring terms of order k/V, or smaller, in equation ( |111|) , we find 



= v^/ + -^ (115) 

Vy/25 



Combining these we find that a + b = —l/^/r, or 6 = ^ up to corrections of order k/V. 

Finally, we can derive an equation for A. Varying A and requiring that equation (|105|) 
be stationary gives 

? 1 

2tX = J dypix)^^ (116) 

a 

It is possible in general to evaluate the following integral 

b 

dyp{x)^— (117) 

y-x 

The result is 



(A + a)(A + 6)(2r , ^ - 2r + ^^-( , ^ - ^ )) 

y'(A + a)(A + 6) Vx + X ^(a-x)(6-x) ^(A + a)(A + 6) 

(118) 

Taking a limit as x goes to A in the above equation, and substituting into equation (|116D , 
we find 

(A + a)(A + 6)3/2 ^ ' 

Finally, using the consistency and volume equations, as well as the equation for A, we 
find that b = V^,a = -j^ + ^2^3(^)2/3^ and A = ^ - |-i.2i/3(^)2/3. 



XVIII. APPENDIX B: CONTRIBUTION TO Sk OF ORDER V/VK 

As discussed in section XI, the calculation of the fluctuations in Aj and the calculation 
of the determinant of Mij is only approximate. There will be corrections to this calculation 
which will produce correction to Sk of order V/^/k. 

We will not precisely evaluate these terms, but we will show that they are of order V/y/k; 
this is slightly surprising since there are many other corrections to Sk which are of order V/k. 
For example, the calculation of the integral over q, so that equation (0), which includes 
only Gaussian fluctuations in c, is not quite right. If calculated more carefully, there would 
be corrections of order V/k from the cubic and quartic terms in the action q. 
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We will look carefully only at the corrections to order V/yk^ the correction to the 
calculation of the determinant of Mjj. The calculation will only be sketched, to save space. 

XI, are equal to 



Off-diagonal terms in Mj^ after removing terms in Mjj of order (3"^ as discussed in section 



2E E ^^r^K^l (120) 

H=l a=k+l ^o- 

These terms have mean square \{Gtj) > where G^^ is the Greens function of matrix H^ + X. If 
one writes the determinant as e'^'°§^'-', and perturbatively expands the log in the off-diagonal 



V 



i2 



2 Z^ k\^ijJ 

terms, the first correction one finds is e '•^=^ . This might make one think that these 

corrections will only change Sk by order V/k. However, although in the paramagnetic phase 
(G^) is of order 1/V this is not true in the spin glass phase, and the perturbative expansion 
diverges and is therefore not valid. 

We can use a trick to get around this. The desired determinant can be written 

'^•^ (121) 



(det(i7> + A))'= 

Here, /x = l...k and a = k + 1...V. The vectors v"' are the eigenvectors of (iJ> + A). In the 
basis of these eigenvectors, (iJ> + A) is diagonal. For simplicity, let us evaluate this with a 
specific choice of 0f . Let us take (p} = P ior i = l...V/k and cf)} = ioi i = 1 + V/k...V. 
Also, we will take (j)j = j3 for i = 1 + V/k...2V/k and zero for all other i. We will follow 
this pattern for all other fi, so that for a given i, 0f is non-zero only for one value of /i. 
The matrix i7> + A has V ~ k non-zero eigenvalues, and can be thought of as a (V — k)-hy- 
{V — k) matrix. Using the particular choice of 0f here, the numerator of the above equation 
can be written as the product of k determinants, one for each /i. Each determinant is the 
determinant oi a. {V — k — V"/A;)-by-(V^ — k — V/k) random matrix since for the given choice 
of 0, for each fj, the factor HilZ] 2%l)^v'^(j)':ll(j)'^v'^ip'/^) removes V/k different factors of ip^^ip^^ 

IJ,,U 

from the integral. In fact, this result is not specific to this choice of 0f , but is easiest to 
obtain with this particular choice. 

Now, if we evaluate the product of the determinants of the k different V — k — V/k-hj- 
V — k — V/k random matrices, and divide by the k-th power of the determinant of iJ> + A, 
we find that the result is {2^/T)^ , times some corrections which change Sk by terms of order 
V/vk. The calculation of the ratio of these determinants is elementary and will not be done 
here. It is simply a matter of using the Wigner semicircle distribution for the eigenvalues of 
the different matrices, and evaluating the determinant by taking a product of eigenvalues. 

Similarly, there are corrections to the fiuctuations in Aj, involving the contribution to 
these fiuctuations resulting from {det{H^ + Aj))^. Looking at a perturbative expansion of 
the determinant similar to that performed above, one would again think that the corrections 
to the fiuctuations change Sk by order V/k, but again the perturbation expansion involves 
{GijY and is invalidated in the spin glass phase. One again finds corrections to Sk of order 
V/Vk. 
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Finally, the corrections to Sk discussed here depend on the average value of Aj. This can 
lead to a shift in A away from the value calculated in appendix A. These corrections can 
adjust the value of the gap, A + a, but for k oc V"^^^ they will not alter the scaling of the 
gap with k; if we calculate an effective action for A combining the corrections discussed here 
with the results from appendix A, all the terms will be of the same magnitude and X + a 
will still scale as (y)"^^^- 
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FIG. 1. Illustration of the construction of a series of systems saturating the bound. 




FIG. 2. Wigner semicircle distribution of eigenvalues. The horizontal axis is energy, the vertical 
axis is eigenvalue density. The semicircle extends from 



^ to ./i. 




FIG. 3. Altered level distribution of H + X in large A^ system. The spike contains k eigenval- 
ues, the smaller semicircle contains V — k eigenvalues. Superimposed is a semicircle of a matrix 
containing V eigenvalues. The horizontal axis is energy, although the scale is shifted so that the 
spike is at zero energy. 
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